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Using an effective Hamiltonian including the Zeeman and internal interactions, we describe the 
quantum theory of magnetization dynamics when the spin system evolves non-adiabatically and 
out of equilibrium. The Lewis- Riesenfeld dynamical invariant method is employed along with the 
Liouville-von Neumann equation for the density matrix. We derive a dynamical equation for mag- 
netization defined with respect to the density operator with a general form of damping that involves 
the non-equilibrium contribution in addition to the Landau-Lifshitz- Gilbert equation. Two special 
cases of the radiation-spin interaction and the spin-spin exchange interaction are considered. For 
the radiation-spin interaction, the damping term is shown to be of the Gilbert type, while in the 
spin-spin exchange interaction case, the results depend on a coupled chain of correlation functions. 

PACS numbers: 76.20.-(-q, 72.25.Ba 



I. INTRODUCTION 

Magnetization dynamics in nanomagnets and thin 
films is rich in content, including such phenomena as 
giant magnetoresistance spin-current-induced mag- 
netization reversal and adiabatic spin pumping 
The study of the magnetization dynamics is motivated 
by theoretical interest in a deep and fundamental under- 
standing of the physics of magnetic systems on a short 
time scale and out of equilibrium. 

The results are important in gaining an understand- 
ing of the technological applications of magnetic systems 
to areas such as high-density memory and data storage 
devices. Recent results, in real time, provide the behav- 
ior of spins in particular set-ups of magnetic fields and 
measuring the spin flips as a function of time 01 j @ ■ 

The Landau-Lifshitz-Gilbert (LLG) equation 
provides a plausible phenomenological model for many 
experimental results. Recently, the LLG equation and 
the Gilbert damping term have been derived from an 
effective Hamiltonian including the radiation-spin inter- 
action (RSI) 0] . It has been assumed there that the spin 
system maintains quasiadiabatic evolution. 

However, a magnetic system whose Hamiltonian Ti.{t) 
evolves nonadiabatically, i.e. in a non-equilibrium state, 
deviates far from quasi-equilibrium, and its density op- 
erator satisfies the quantum Liouville - von Neumann 
(LvN) equation 



ih 



dp 
at 



0. 



(1) 



In the present work, we aim to derive a magnetization 
equation including a damping term for such a nonequilib- 
rium magnetic system. To find its nonadiabatic quantum 
states, we employ the Lewis-Riesenfeld (LR) dynamical 
invariant method (lo| . This method originally designed 
for the nonequilibrium evolution of time-dependent quan- 
tum systems has been successfully applied to a variety 
of problems, including the nonadiabatic generalization 
of the Berry phase for the s pin dyiiamics, nonequilib- 
rium fermion systems etc 11 1, 12|,lij]- l3> the time- 
dependent invariants have been used for constructing the 
density operator for nonequilibrium systems. Following 
[1^ , we construct the density operator and then use it to 
determine the time evolution of magnetization. 

The damping term represents magnetization relaxation 
processes due to the dissipation of magnetic energy. Var- 
ious kinds of relaxation processes are usually melded to- 
gether into a single damping term. Relativistic relaxation 
processes result in the Gilbert damping term with one 
damping parameter, while for the case of both exchange 
and relativistic relaxation the darnping term is a tensor 
with several damping parameters [l5|, [l6[. 

The relaxation processes are specified by interactions 
of spins with each other and with other constituents of 
the magnetic system. A derivation of the damping term 
from first principles should therefore start with a micro- 
scopic description of the interactions. Even though such 
microscopic derivation of damping has been performed 
for some relaxation processes (for instance, [131 a full 
version of derivation for the Gilbert damping term has 
not yet been given, and in particular for a system in 
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non-equilibrium. 

Herein, we consider a general spin system without 
specifying the interaction Hamiltonian and related relax- 
ation processes. We start with a system of spins precess- 
ing in the effective magnetic field He// neglecting for a 
moment mutual interactions. Then at a fixed time later 
interactions in the system are switched on and influence 
the original precessional motion. 

The interactions are assumed to be time-dependent, 
and the spin system evolves nonadiabatically out of equi- 
librium trying to relax to a new equilibrium magnetiza- 
tion. We perform a transformation, which is analogous to 
the one used in the transition to the interaction picture, 
to connect the density and magnetic moment operators 
before and after the time, when a new nonequilibrium dy- 
namics starts, and to find an explicit expression for the 
interaction contribution to the magnetization equation. 

Our paper is organized as follows. In Sec. II, the mag- 
netization equation for the system of spins with a gen- 
eral form of interactions is derived. In Sec. Ill, the case 
of the radiation-spin interaction is considered, which is 
shown to produce the Gilbert damping term even for sys- 
tems that are not in equilibrium. The contribution of the 
non-dissipative part of the radiation field to the magneti- 
zation equation and magnetization algebra is discussed. 
Sec. IV focuses on a special type of spin-spin interac- 
tions. We conclude with discussions in Sec. V. 



II. MAGNETIZATION EQUATION 

Let us consider a quantum spin system defined by 

n^Ho + XHi, (2) 

where Hq is the Zeeman Hamiltonian describing the in- 
teraction of spins with an effective magnetic field 



We introduce next the magnetic moment operator 



7io = -7E^'-He//W, 



(3) 



7 being the gyromagnetic ratio and Si being the spin 
operator of the ith atom, while the Hamiltonian 7i/ de- 
scribes the internal interactions between the constituents 
of the spin system, including, for instance, the exchange 
and dipolar interactions between the atomic spins, as well 
as higher-order spin-spin interactions. 

In Eq.(l3]), the effective magnetic field Heff is given 
by the energy variational with magnetization , VLeff — 
—5E{M) j 5M, where E{M) is the free energy of the mag- 
netic system. This field includes the exchange field, the 
anisotropy field, and the demagnetizing field, as well as 
the external field, tiext- 

The interaction terms included in Tij are in general 
time-dependent, being switched on adiabatically or in- 
stantly at a fixed time to- The parameter A in (21) can 
be chosen small in order to take into account the higher 
order effects perturbatively. 



M 



sn 



(4) 



which is the response of the spin system to the exter- 
nal field. The magnetization is defined as an ensemble 
average of the response 



1 



M^{M) = -Tr{pA^}, 



(5) 



where p is the density operator satisfying the LvN equa- 
tion ll]) and V is the volume of the system. The explicit 
form of the density operator will be shown below. 

For systems in equilibrium, the Hamiltonian itself sat- 
isfies the LvN equation and the density operator is ex- 
pressed in terms of the Hamiltonian. For nonequilibrium 
systems, whose Hamiltonians are explicitly time depen- 
dent, the density operator is constructed by making use 
of the time-dependent adiabatic invariants. 



A. Zeeman precession 

Let us first derive the magnetization equation for 
the Zeeman Hamiltonian. If the interactions, 7i/, are 
switched off, the magnetic moment operator and the 
magnetization are 



and 



Mo = (>Jo>o = ^^Tr{po5.} 



(6) 



(7) 



respectively, the operators Mq, a = 1,2,3, fulfilling the 
SU{2) algebra 



(8) 



where the summation over repeated indices is assumed. 
In Eq.©, the subscript "0" in the symbol (...)q indicates 
using of the density operator po, which satisfies the equa- 
tion 



ot 



(9) 



Let T()(t) be a non-trivial Hermitian operator, which 
is a dynamical invariant. That is, To{t) satisfies the LvN 
equation 



dio dio 1 



dt 



dt 



2^0, ^o]_ = 0. 



(10) 



As shown in [T3| , the eigenstates of Xq {t) can be used for 
evaluating the exact quantum states that are solutions 
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of the Schrodinger equation. The linearity of the LvN 
equation allows us to state that any analytic functional 
of Xo{t) satisfies the LvN equation provided that To{t) 
satisfies the same equation. In particular, we can use 
io{t) to define the density operator po as [Ijl 

Mt) = ^ e-^^'>^'\ Zo - Tr{e-'5*°«}. (11) 

Here /3 is a free parameter and will be identified with the 
inverse temperature for the equilibrium system. 
The LvN equation for Tq is formally solved by 



Io{t) = U{t,to)Io{to)U{to,t), 



(12) 



where 



J7(io,t) = Texp{^^ rfT7io(T)}, (13) 

and T denotes the time-ordering operator. 

As the Zeeman Hamiltonian is linear in spin operators, 
we can take lo (t) of the same form |18| 



(14) 



where Rq is a vector parameter to be determined by a 
dynamical equation. Then, Eq.® becomes 



dt 



-7RoxHe// =0 (15) 



and we obtain an equation for the vector parameter given 

by 



dRo 



-|7|Ro X He//. 



(16) 



This equation explicitly describes the vector precessing 
with respect to the field He//. Thus, without loss of 
generality, we can identify Rq with the magnetization 
vector Mq (up to a dimensional constant factor) and 
Eq. pop with the equation of motion of magnetization 
for the Zeeman Hamiltonian Ti.Q. 

An alternative method to obtain the magnetization 
equation is just to differentiate both sides of Eq.(I7]) with 
respect to time and use Eq.([H]). In this way, we arrive at 
Eq.ldll) with Ro = Mo- 



B. Interactions and damping 

In the case, when the interactions are present, the den- 
sity operator for the full Hamiltonian may be written as 



p(t) = le-'5*W, Z = Tr{e-''*W}, (17) 



where X is the invariant operator for the system. As 
the interaction Hamiltonian is in general non-linear in 



spin operators, T{t) cannot be taken in the form given 
by Ea. p4)) . Moreover, the form oiT{t) can not be deter- 
mined without specifying the interactions. 

To derive the magnetization equation in this case, we 
proceed as follows. We perform, on p{t), the transforma- 
tion defined by the operator p3)) . 



(18) 



P Pint = U{to, t)p{t)U{t, to), 



removing the Zeeman interaction. For systems with Tig 
constant in time, the operator U{tQ, t) — exjp{(i / h)Ti.Q{t— 
to)} leads to the interaction picture, which proves to be 
very useful for all forms of interactions since it distin- 
guishes among the interaction times. For our system with 
both Tio and Tij dependent on time, the operator (jl3[) 
plays the same role, removing the unperturbed part of 
the Hamiltonian from the LvN equation. 
Substituting Eq.lUHl) into yields 



ih 



dpi: 



dt 



X 



'Hint, Pint 



where 



n,nt{t) = u{to,t)ni(t)u{t,to) 



(19) 



(20) 



The magnetic moment operator and the magnetization 
become 



where 



and 



Mi = -A 



M 



<5He 



Tr{p„f(A^O,jnt +Ml^int)}, 



(21) 



(22) 



(23) 



Mo.int and Mi^int being related with and Mi, re- 
spectively, in the same way as Tiint is related with Tii. 

The operators A^o.int, Mi^int are generally used to 
calculate the magnetic susceptibility [l^. Let us show 
now how these operators determine the time evolution of 
magnetization. The evolution in time of A^o, int is given 
by the equation 



dMo, 



dt 



-U{tQ,t) 



no, Mo 



U{t,to) 



iMo^int X He//, 



(24) 



which is analogous to Eg. p^ . It describes the magneti- 
zation precessional motion with respect to He//. 
The equation for Misnt, 



dM 



I .int 



dt 



h 



U{to,t) 



Ho, Ml 



U{t,to) 



4 



(25) 



describes more complex magnetization dynamics gov- 
erned by the interaction Hamiltonian Tij. In addition, 
J^i depends on time explicitly. However, this dynamics 
includes the precessional motion as well, sine the inter- 
actions induced magnetization is a part of the precessing 
total magnetization M. Introducing 



Ho, Mi 



^Ml X He// 



(26) 



to represent deviations from the purely precessional mo- 
tion, we bring the equation for Mi.int into the following 
form 



dt 



= iMl^int X He// 



+ ?7(to,i)(^+D/)j/(i,to). 

Taking the time-derivative of M given by Eq. 
using Eq. p9|) . we obtain 



(27) 
and 



dm 

Hi 



V 

x_ 

ih 



^Tr{p„,( 



dM dMi ,int 



dt 



dt 



■MoAnt + -Ml Ant, 'Hint 



(28) 



Substituting next Eqs.jSl]) and ^ into Eq.(l28]), finally 
yields 



^ = -|7|MxHe//. 



D, 



where 



D 



M,ni 



,dMi 
' dt 



(29) 



(30) 



Therefore, Eq. p9|) is the magnetization equation for the 
system specified by ^ . This equation is general since it 
is derived without specifying TLi. The D-term contains 
all effects that the interactions, TtI/, can have on the mag- 
netization precession, so that Eg. ipS]) is complete. 

The contribution of TLi to the D-term in the magne- 
tization equation can be divided into two parts. One is 
proportional to ([A^,7i/]_) and is related to the change 

in the density matrix when the interactions of 7i/ are 

switched on. The second part (^^^ + D/) originates 
from the change in the magnetization itself. Which part 
of D is dominating depends on the nature of the interac- 
tions. 



III. RADIATION-SPIN INTERACTION 

One of the important issues in the study of the mag- 
netization dynamics is the relaxation phenomena. The 
magnetization relaxation mechanism can be introduced 
by various interactions such as spin-orbit coupling and 
two-magnon scattering. In this section, we calculate the 
D-term for the magnetization relaxation process, which 
is induced by the RSI. 



A. Dissipative radiation field 

The RSI approach is an effective field method, in which 
contributions to the magnetization relaxation are effec- 
tively represented by the radiation-spin interaction [oj. 
In this method, the damping imposed on the processing 
magnetization originates from the magnetization preces- 
sional motion itself. It is assumed that a radiation field 
is induced by the precession and that this field acts back 
on the magnetization producing a dissipative torque. 

The Hamiltonian for the RSI is 



Hf 



where 



H_ 



A(M X He// - aM^He//) 



(31) 



(32) 



is the dissipative part of the effective radiation field, i.e. 
the part responsible for the dissipative torque, and M is 
the magnitude of magnetization. The parameter a will 
be specified below, while A in Eq. ([2]) is now the radiation 
parameter. 

The RSI contribution to the magnetic moment is 



Mia = -KaM^Mo-Mo x M), 
and the total magnetization becomes 



(33) 



M = {l-XaM^){Mo)+X{MQ) xM. (34) 

It turns out that the vectors M and (Ma) arc parallel. 
Indeed, taking the vector product of both sides of Eq, 
with (A^o): yields the equation 

(Mo) X M = X{Mq) X {Mq) X M 



= X{Ma) (l- XaM^){Mo) -M , (35) 
which is valid only if M is parallel to (Mq), so that 

A straightforward calculation shows that 

1 r 



{Di) = -A( 



ih 



M,ni 



5 



7A 



1 - XaAP 



M X M X H 



(37) 



so that the contributions of D/ and [A4, Hi]_ to Eg. ([30]) 
cancel each other. Calculating next the time derivative 
of A^/,d and using Eq. p6|) . we obtain 



D = aM X 



dM 

IT' 



(38) 



provided the following relation between the parameters 
a and A holds 



A 



1 - AaM2 ■ 



(39) 



Therefore, Eq. p9|) takes the form of the LLG equation, 
a becoming a damping parameter. 

The relation given by Eq. ([55]) reflects the origin of the 
radiation field. Both the radiation and damping parame- 
ters depend on the magnetization. The dimensionless pa- 
rameters independent of M are Ao = AM and ao = aM 
with the relation 



Ao 



1 + al 



(40) 



The equivalent form of the LLG equation, which is 
more suitable for calculations, is 

^ = -|7|M X [(1 - XaM^)ll,ff + AM X H,;/] . (41) 

Let us assume that tlef / is a uniform static field in the z- 
direction, i.e. He// = (0,0, iJz). Then the magnetization 
equation becomes, in component form. 



^M, = Awo Mp, 
at 



(42) 
(43) 



where = |7|^^z is the frequency of magnetization pre- 
cession and = M^+M^ = A'P-M^. These equations 
are solved exactly by 



and 



M, = M ■ 



Mr, 



tanh[AocJo(i — ^o)] + d 
1 + d ■ tanh[Aoti^o(^ ~ ^o)] 



(44) 



MVl - d"^ 
cosh[AoWo(^ - to)] 



(45) 



1 + d ■ tanh[Aoa;o(t — to)] ' 

where d stands for the initial condition 

A4(t^M 
M ■ 

In the limit t — *■ oo, Mz tends to M and Mp tends to zero, 
so that during the relaxation process the magnetization 
vector tends to be parallel to the effective magnetic field, 
and the relaxation characteristic time r is l/(AoCi-'o)- 



B. General radiation field 

The RSI approach can be directly used for studying the 
effect of the real radiation-spin interaction in the mag- 
netization relaxation process. In that case, one has to 
distinguish the real radiation field contribution to the 
damping parameter from the effective one. It is there- 
fore important to consider possible generalizations of the 
ansatz given by Eq. ([5^ to get a more detailed picture of 
the RSI. 

The dissipative radiation field in Eq.(l32|) is a compos- 
ite of two fields, which are parallel and perpendicular to 
He//, respectively. The field parallel to He// changes 
only the frequency of the magnetization precessional mo- 
tion, while the torque He// x M introduces a damping 
effect as well. 

The choice of the field llf made in Eq. (|5^ is not a 
general one. Let us modify it by applying additional 
fields and see how the parameters in the magnetization 
equation are changed. If we apply an additional field in 
the direction of He/ / and modify Hf as 



HJ! 



H„ — AnH 



one//, 



(46) 



where Ao is an arbitrary dimensionless parameter, then 
this results in a change of the magnitude of the effec- 
tive magnetic field and therefore in the change of the 
frequency of the magnetization precession as 



LUo 



uj = wo(l - Ao). 



(47) 



For Aq > 1, the magnetization vector turns upside down 
and the precession continues in an opposite direction. 

The magnetization damping is not affected by the ad- 
ditional field. Rescaling of He// in Hf leads to a change 
of the radiation parameter as well. 



Ao 



1 



1-Ao 



=^Ao, 



(48) 



so that the relaxation characteristic time remains the 
same. 

If we apply an additional torque, 

Hf ^Hf + AoMxHe//, (49) 
then the parameters in the ansatz ([5^ change as follows 
Ao ^ Ao + Ao, (50) 



ao 



Ao + Ao, 
Ao 



-ao, 



Ao + Ao 

and the relaxation characteristic time becomes 



T 1 



Ao^-^ 



Ar 



(51) 



(52) 



decreasing for Ao > and increasing for — Ao < Aq < 0. 
If the additional torque is stronger than the original one 
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FIG. 1: The dependence of the z-component of the magne- 
tization vector, Mz/M, on time x = {t — to)/r is shown for 
different values of a and for d — 0.4. 



1.2 

1 

0.8 
0.6- 
0.4 
0.2- 











V ^ 


a=0 




- - - .3=0.5 






— - ■a=-0.5 






a=-1.5 






• \ \ \ 









12 



FIG. 2: The time dependence of the perpendicular component 
of the magnetization vector, Mp/M, is plotted for different 
values of a and for d = 0.4. 



and in the opposite direction, i.e. Xq < — Ag, then the 
magnetization vector turns over again. For Aq — — Aq, 
the additional and original torques cancel each other, and 
there is no damping. 

In Fig[T]and Figl2]we plot the solutions for and Mp 
given by Eq. (|44|) and Eq.((45|l. respectively, in the pres- 
ence of an additional torque (nonzero a's) and without it 
(a — 0), where a = Xq/Xq. We observe a change in the 
relaxation time depending on the sign of a. For a > 0, 
the additional torque is in the same direction as the orig- 
inal one, while for a < their directions are opposite. 
A drastic change occurs for a < — 1, when tends to 
{—M) indicating the overturning of M. The magnetiza- 
tion vector first becomes perpendicular to He//, when 
Mp/M reaches its maximum value, and then it turns up- 
side down. 

The radiation field can contain a non-dissipative 
part as well, i.e. H^. = HJ?-f H". Let us assume that H" 



is parallel to M and we take it of the form 



H;^ = «A^(M.He//), 



(53) 



where k is an arbitrary parameter. 

The field given by Eg. ([55)1 does not change the mag- 
netization precessional motion. However, the magnetic 
moment operator is given as 

M^M+Mi,n=M + KX^iMo-M). (54) 

The relation between the vectors (A^o) and M becomes 

M 



{Ma 



1 - X{aM + k)M' 



(55) 



Proceeding in the same way as before, when the non- 
dissipative part of the radiation field was omitted, to cal- 
culate the total radiation field contribution to the D- 
term, we obtain the same Gilbert-type structure of the 
damping term with the same damping parameter, i.e. the 
part of the radiation field that is parallel to M does not 
contribute to the magnetization equation. 



C. Magnetization algebra 

The magnetization algebra, i.e. the algebra of mag- 
netic moment operators, gives a further insight into the 
RSI. Let us study how the non-dissipative part of the 
radiation field contributes to this algebra. 

Although we can omit the non-dissipative radiation 
field in the procedure obtaining the equation of motion 
for magnetization, it is important to take this field into 
account when we construct the magnetization algebra. 
Indeed, the operators M"" = Mq + Mj,i without the 
non-dissipative radiation field contribution fulfil the al- 
gebra 



ih"fe'""'{{l - XaM') ■ M' 



-XM ■ MjJ, 



(56) 



which, however, contains Mi^n on its right-hand side and 
therefore the algebra is not closed. 

Let us define, as in Eq. (f54|) . the total magnetic moment 
operator Aitot = M + A^/,,i. Then the magnetization 
algebra becomes 



Ml,, Ml, 



ih-fe'''"'{{l ~ XaM^ + kXM) ■ 



+ {-XM - 1 + XaM^) ■ MIJ. (57) 
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Fixing next the parameter n by taking it as a solution of 
the equation 



1 

K 

K 



2 



that is 



K± 



- — (lTx/l + «g), 
ao 



(58) 



(59) 



brings the algebra into the closed, standard form for 
SU{2) algebra 



injtote''''^M^tot 



(60) 



The non- vanishing commutator [A4o, Hi] _ in Eq. pOl) 
is the only contribution of the spin-spin interaction to 
the magnetization equation, resulting in 



where 



and 



D 



A, 



-M X J7, 



Ml,Ml =MIMI + MIMI 



(64) 



(65) 



(66) 



where 



Itot = ±7 



(61) 



For < ao < 1, > and k_ < 0, the sign (+) in 
Ea. (|6ip corresponding to the and sign (— ) 

to K = K_. 

Therefore, the RSI preserves the form of the magneti- 
zation algebra (c.f. Eq.®), by renormalizing the gyro- 
magnetic ratio, only if the interaction of spins with the 
non-dissipative part of the radiation field is included in a 
proper way. If the direction of H" is the same as that of 
M and n — k^, the gyromagnetic ratio remains negative, 
while its absolute value decreases. If the direction of H" 
is opposite to that of M and k = k_, then the gyromag- 
netic ratio changes its sign and becomes positive. 



IV. SPIN-SPIN INTERACTIONS 

The spin-spin interactions among the spins in the sys- 
tem introduce many body effects, which can be treated 
perturbatively in the weak coupling regime. In this case 
the D-term can be expanded in powers of A. To demon- 
strate this, we consider the spin-spin interactions of a 
specific type. The interaction between spins is usually 
an exchange interaction of the form 



2 J ^ ^ SiS j 

i,3 



—Mo, 

7 



(62) 



the coupling constant J being called the exchange inte- 
gral. We generalize the ansatz given by Ea. ([62)l by as- 
suming that the exchange integral depends on the mag- 
netization and introduce the spin-spin interactions as fol- 
lows 



(63) 



where J"'' = XM"-M^. Since fii does not depend explic- 
itly on the external field, its contribution to the magnetic 
moment operator vanishes, Mj = 0. 



The correlation function G"'' 
of spin correlation functions. 



is the sum 



(67) 



excluding the self-interaction of spins. For the standard 
ansatz given in Eg. ([5^ . D = and the magnetization 
equation does not change. 

If the spin-spin interactions are turned on at t = tg, 
so that /5(<o) = Po{to), then, integrating both sides of 
Eq.dinD, we find 



A f 

P\{t) = Po(io) + -7 dr 



'H\{t),px{t) 



(68) 



Substituting Eg. ((68)) into the definition of G"'', yields the 
equation 

G''\t) = Gg''(tn) 



- / dTj'''^{T)(e'"'''G'''"^{T)+e'"'''G'""^{T)\, (69) 

7 t/io 



where 



Gah 
= 



/O' 



(70) 



which relates G""^ to the third order correlation function, 
i.e. the correlation function of the product of three mag- 
netic moment operators. 



G 



abc 



(71) 



The correlation function G"*"^, in turn, is related to the 
fourth order correlation function and etc., and we have 
therefore an infinite number of coupled equations for spin 
correlation functions. For any practical calculation this 
infinite hierarchy has to be truncated. That then defines 
the approximation scheme which may be considered on 
the basis of the physical requirements for the system. 
The approximation scheme will depend on the physical 
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properties such as density and on the strength of the 
interactions. 

If the Hamiltonian 7i/ is a smah perturbation to the 
original Hq, wc can solve Eq. (j68p perturbatively. In the 
lowest, zeroth order in A, we replace px(t) by po{ta), so 
that G"'' « Gfito). We choose the initial value for G"'' 
as 

Y^Y.i^fS^)^^!-" (72) 

with /^^ = = 0, /^^ = / and I"'' = for a 7^ 6. We 
also define again the z-direction as the direction of the 
effective magnetic field that is chosen uniform and static. 
Then the D-term becomes, in component form, 

I?, = 2XjIMyM,, (73) 
Dy = -2XjIM,M„ (74) 
= 0. (75) 

producing two effects in the magnetization equation: 
the magnetization is now precessing with respect to 
{Hz — 2X1 Mz), its z-component remaining constant in 
time, {d/dt)Mz ~ 0, and the frequency of the precession 
is 

tUo = wo(l - A^^-^y (76) 

Therefore, in the lowest order of perturbations, when the 
D is linear in A, it shifts the direction and the frequency 
of the precessional motion without introducing damping 
effects. To find a role of the higher powers of A in D and 
to determine how they affect the magnetization equation, 
a truncation of the chain of spin correlation equations is 
needed. This would require a consistent perturbation 
approach to the hierarchy of the coupled equations for 
the correlation functions. 



each other and with other degrees of freedom of the sys- 
tem is specified. For the interactions related to the re- 
laxation processes, this equation provides a general form 
of magnetization damping. 

This paper uses the dynamical invariant method in- 
troduced by Lewis and Riesenfeld It extends its 
applicability to magnetic systems that arc in a non- 
equilibrium state i.e. the various components in the 
defining Hamiltonian are time-dependent. This is in con- 
trast to the earlier attempts to sove a problem for quasi- 
adiabatic evolution. 

The D-term in the magnetization equation has been 
obtained without using any approximation scheme. It is 
exact, accumulates all effects of the internal interactions 
on the magnetization precessional motion and can be a 
starting point for practical calculations. We have eval- 
uated the D-term in two special cases, the RSI and the 
spin-spin interactions. For the RSI, which is linear in 
spin operators, it takes the form of the Gilbert damping 
term, the damping and radiation parameters being inter- 
related. For the spin-spin interactions, it is determined 
by the spin correlation functions, which fulfil an infinite 
chain of equations. A further analysis of the D-term re- 
quires an approximation scheme to truncate the chain in 
a consistent approach to higher order calculations. 

In our work, we have considered a specific type of the 
spin-spin interactions, which do not contribute to the 
magnetization algebra. However, if spin-spin interactions 
depend explicitly on the external field, the form of the 
algebra can change. In this case, the total magnetic mo- 
ment operator becomes nonlinear in M.Oi and this re- 
sults in the magnetization algebra with an infinite chain 
of commutation relations. The chain has to be truncated 
in a way consistent with the truncation of the chain of 
equations for the spin correlation functions in the same 
approximation scheme. 



V. DISCUSSIONS 

We have derived a general form of magnetization equa- 
tion for a system of spins precessing in an effective mag- 
netic field without specifying the internal interactions. It 
can be applied in the study of magnetization dynamics 
of any type, including nonequilibrium and nonlinear ef- 
fects, provided the interaction of individual spins with 
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